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ABSTRACT 

We propose the Apparent Shrinking Criterion (ASC) to interpret the spatial extent, 
Rw, of transmitted flux windows in the absorption spectra of high-z quasars. The 
ASC can discriminate between the two regimes in which corresponds either to the 
physical size, i?Hii, of the quasar H ii region, or to the distance, R™'^^, at which the 
transmitted flux drops to = 0.1 and a Gunn-Peterson (CP) trough appears. In the flrst 
case (H ii region (HR) regime), one can determine the IGM mean H i fraction, xhi; 
in the second (proximity (PR) regime), the value of R^ allows to measure the local 
photoionization rate and the local enhancement of the photoionization rate, Tq, due to 
near by /intervening galaxies. The ASC has been tested against radiative transfer+SPH 
(Smoothed Particle Hydrodynamics) numerical simulations, and applied to a sample 
of 15 high-z {z > 5.8) quasar spectra. All sample quasars are found to be in the 
PR regime; hence, their observed spectral properties (inner flux profile, extent of 
transmission window) cannot reliably constrain the value of xhi ■ Four sample quasars 
show evidence for a local enhancement (up to 50%) in the local photoionization rate 
possibly produced by a galaxy overdensity. We discuss the possible interpretations and 
uncertainties of this result. 

Key words: cosmology: theory - radiative transfer - methods: numerical - intergalac- 
tic medium - cosmology: large scale structure of Universe - quasars: general 



1 INTRODUCTION 

All high-redshift quasar spectra with a complete GP trough 
(Gunn & Peterson 1965) exhibit transmitted flux in an ex- 
tended region between the quasar redshift and the red side 
of the GP trough (Fan et al. 2006, F06). These transmission 
windows provide a powerful, and so far unique, probe of the 
physical state of the InterGalactic Medium (IGM) beyond 
redshift z > 5.7, where the onset of the GP trough erases 
the information encoded in the Lya forest. 

The detection of transmitted flux at such high redshift 
provides a solid evidence for a reduced neutral hydrogen 
fraction (2;hi) within distances of several physical Mpc0from 
the quasar (F06). Such H i decrease can be attributed to 
the locally enhanced photoionization rate produced by the 
quasar ultraviolet radiation. However, the physical interpre- 
tation of such transmission windows in terms of the IGM 
properties is not straightforward. Are these transmission 
windows the spectral counterparts of quasar H ii regions 
(HR), with the red side of the GP trough corresponding 



^ Unless otherwise stated, distances are given in physical units. 



to the location of the quasar ionization front? Or are they 
instead proximity regions (PR) resulting from the reduced 
optical depth in the vicinity of a quasar whose H ii region 
extends further into the IGM ? We will refer to these two in- 
terpretations as the HR and PR regime, respectively. It has 
been demonstrated that both regimes can be used to repro- 
duce the observed spectra, depending on the assumed quasar 
luminosity and lifetime, and the ionization and thermal state 
of the surrounding IGM (Bolton & Haehnelt 2007a; Maselli 
et al. 2007 [M07]; Wyithe, Bolton & Haehnelt 2008 [W08]). 

Previous studies have shown that under the assumption 
of the HR regime, the relatively small extent of the observed 
transmission windows detected in current data would imply 
a significant neutral hydrogen fraction still present at z ~ 6, 
, that is xni > 0.1 (Wyithe & Loeb 2004; Wyithe, Loeb & 
Carilli 2005; F06). These results are supported by indepen- 
dent analysis which do not explicitly assume the HR regime, 
based either on fitting the probability distribution function 
(PDF) of the Lya optical depth, r, within the transmis- 
sion windows (Mesinger & Haiman 2007), or on the com- 
parison of Lya and Ly/3 absorption spectra (Mesinger & 
Haiman 2004). However, it is difficult to reconcile these re- 
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suits with the analysis of dark gaps in quasars/gamma-ray 
bursts absorption spectra (Gallerani et al. 2008a, 2008b), 
as well as with other independent probes of the hydrogen 
reionization history, as the Lyman Alpha Emitters surveys 
[e.g. Malhotra & Rhoads 2006; Dawson et al. 2007), and 
the integrated electron scattering optical depth measured 
by WMAP5 data release (Komatsu et al. 2009). All these 
studies suggest a;Hi<S 1 at z « 6 and are consistent with the 
detection of transmitted flux in high-redshift quasar spec- 
tra provided that the PR interpretation of the transmission 
windows applies. In this regime, the spectral extent of the 
windows is set by the onset of the GP trough within the 
quasar H ii region. In M07 this effect has been termed ap- 
parent shrinking, to indicate the possibility that the physical 
extent of the transmitting window, Rw , does not correspond 
to the actual H ii region size, -Rhii, and in fact underesti- 
mates it. 

The condition < Rmi is likely to arise from the 
combination of two factors: (i) the average IGM density gets 
high enough to produce complete Lya absorption in spite 
of the high gas ionization level, and (ii) the quasar photo- 
ionization rate drops below the minimum value required to 
preserve IGM transparency inside the quasar H ii region. 

To discriminate between the HR and PR regimes we 
propose here a simple (but yet solid) criterion, the Apparent 
Shrinking Criterion (ASC). Such method, in addition, allows 
to determine xm{z) in the HR regime, or to quantify local 
deviations from the mean cosmic photoionization rate along 
the line of sight (LOS) to the quasar if the PR regime holds. 



2 THE APPARENT SHRINKING CRITERION 

In this Section we introduce the basic idea underlying the 
Apparent Shrinking Criterion, whose main use is to discrim- 
inate between the HR and PR scenarios described above by 
using a given high-redshift quasar spectrum. The ASC is 
based on the comparison among three spatial scales: (a) the 
physical size of the spectral transmission window, R^, (b) 
the size of the quasar H ii region, -Rhii, and (c) the max- 
imum value that R^ can take, TC^""^ . These are discussed 
separately in the following. 



2.1 Size of the transmission window 

R^ is a directly measurable quantity in high-z quasar spec- 
tra and it denotes the physical extent of the spectral region 
close to the quasar in which a non-zero flux can be detected. 
The determination of this quantity is affected by uncertain- 
ties on quasar redshift, zq, and on the flux detection limit of 
the instrument. An additional complication is represented by 
the fact that the IGM transmissivity at 2 ~ 6 is a mixture of 
dark gaps and transmission peaks, which makes the identifi- 
cation of the red side boundary of the GP trough uncertain 
(e.g. isolated peaks in the translucent IGM can merge with 
the transmission window region associated to the quasar). 
Here we adopt the operational definition given in F06, which 
states that 7?™ represents the physical distance from the 
quEisar at which, when smoothed by a top-hat filter of 20 A, 
the observed flux drops below a transmission threshold of 
0.1 (corresponding to a Lya optical depth r = 2.3). 



2.2 Size of the quasar H 11 region 

-Rhii denotes the mean radius of the quasar H 11 region. Due 
to the gas inhomogeneities and to radiative transfer effects, 
quasar H 11 regions deviate significantly (up to 15%, e.g. 
M07) from spherical symmetry with a resulting variation 
of the H II region radius along different LOS (see for exam- 
ple Fig. [TJ . Previous studies including RT calculations have 
shown that the radius of a quasar H 11 region expanding in a 
homogeneous IGM approximates well the actual mean value 
of -Rhii along different lines of sight in the inhomogeneous 
case. Hence, we use the following expression to determine 
the average HII region size: 

««-^-KTofeTTfc)"(iii)"^„,"'Mpe, (1, 

The above expression is valid under the additional assump- 
tion that the quasar lifetime, Iq ~ 10^~* yr, is shorter than 
the gas (volume-averaged) recombination timescale, 

t,ec(HI) = 2.6 C-^ ('^) ' Gyr; (2) 

where the clumping factor C = (n^)/(n)^ > 1 is meant to 
include the effects of density inhomogeneities inside the ion- 
ized region; we have used a helium to hydrogen number ratio 
y = 0.08. Due to its low density, typical IGM recombination 
times are long, and exceed tg. 



2.3 Maximum size of the transmission window 

A third fundamental spatial scale entering the formulation 
of the ASC is set by the monotonic dependence of the mean 
transmitted flux, F = e~'^''ff, on the distance from the 
quasar. The detection of the GP trough in all spectra of 
quasars with z > 5.7 implies that at these redshifts the 
mean IGM opacity grows large enough to completely sup- 
press transmission; due to instrumental sensitivity the de- 
tection of GP trough implies an upper limit on the mean 
transmitted flux of -F < 0.1. (see e.g. F06). As a result of 
their large ionizing power, in the vicinity of extremely lu- 
minous quasars the IGM opacity is suppressed with respect 
to average and some flux shortwards of the Lya line can be 
detected. Due to the geometrical dilution of the quasar ion- 
izing flux, -F is a monotonically decreasing function of the 
distance from the quasar, -R, and it will eventually fall below 
the detection threshold -F = 0.1; at even larger distances, the 
quasar photoionization rate drops below the one produced 
by the background sources (QSOs and galaxies). Given the 
quasar ionizing photon rate, N-,, and the mean intensity of 
the ionizing background (or equivalently, a^ni), one can de- 
fine a radius -R™"^, as the distance from the quasar such 
that F{R^'^^) — 0.1, under the hypothesis that the quasar 
fiux extends to infinity. Stated differently, -R™"^ represents 
the maximum spatial extent of the transmission window -R^ 
associated to a given quasar. 

How can we estimate -R™"^? Analogously to -Rhii in 
eq. [1] -R™"^ is calculated in the ideal case of a uniform IGM; 
deviations might occur due to fluctuations along different 
LOS, as well as to uncertainties in the modeling of the IGM 
and quasar properties. These effects will be quantified in the 
next Section. 
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For a uniform IGM the mean optical depth can be ex- 
pressed as a function of redshift and of the neutral hydrogen 
density as follows: 

r = 4.25. 10 h -h:A( — j j j , 

where A = p/{p) denotes the gas overdensity. The mean 
transmitted flux is however sensitive to the probability dis- 
tribution function (PDF) of the density field, -P(A), and can 
be estimated as: 

F = y" P(A)e-"<'^'dA. (4) 

As quasars typically reside in biased (i.e. high-cr peaks) re- 
gions of the cosmic density distribution, -P(A) might differ 
in their vicinity from the one describing the general IGM. 
Under these conditions we do expect in particular that the 
peak of the PDF is shifted towards larger A values. In the 
following, though, we will neglect this bias, as both theoret- 
ical models (Barkana 2002; Wang et al. 2009) and observa- 
tions (Guimaraes et al. 2007) find the ambient overdensity 
to be significant only within a distance (« 1 Mpc) much 
smaller than the observed values of i?^, : the minimum value 
measured to date is the one in the quasar J1623-f3112 at 
z = 6.22 with = 3.6 Mpc. 

In general, F is a function of the distance R from the 
quasar, through the dependence of r on xni, as in eq.|3] The 
neutral fraction, in turn, can be derived from the condition 
of photoionization equilibrium at any given A once the local 
photoionization rate r(_R) is assigned. Close to the quasar, 
r{R) Ri Tq{R) = NjO-h/^tvR'^, with an denoting the 
frequency-averaged photoionization hydrogen cross-section; 
furthermore within distances R ~ -Rhii, r(_R) 2> Fs, the 
mean value in the IGM. The expression for Fq given above 
neglects the flux suppression due to residual H i in the ion- 
ized region; as discussed later on, this is a good approxi- 
mation supported by the results from full 3D cosmological 
radiative transfer simulations. An additional contribution, 
Fg, to the total photoionization rate along the LOS can 
come from galaxies intervening along the line of sight. 

Then the local photoionization rate can be expressed 

as: 

F(i?) =FB+FG(i?) + FQ(i?). (5) 

We initially neglect Fg and assume that FQ(i?) extends to 
infinity (i.e. with no truncation at the H ii region ionization 
front whose location is unknown), and discuss the implica- 
tions of such assumption in Sec. 4. 

In summary, the estimate of the i?™"^ value associ- 
ated to a given quasar requires an a priori knowledge of 
A'^,, Fs and of the mean IGM temperature, T, which en- 
ters the determination of xui (and r, see eq. [Sjl via the 
H II recombination rate. We discuss the choice of these quan- 
tities and the associated uncertainties in Sec. 3. Once these 
quantities are fixed, R^""^ can be estimated from photoion- 
ization equilibrium, and by requiring that F{R^°'^) = 0.1. 
From the analysis of RT calculations, it appears that pho- 
toionization equilibrium is a good approximation in the cases 
of interest. Note that, 7?™"^ is univocally determined by such 
procedure, being t{R) a monotonically increasing function 
of i? as a consequence of the geometrical dilution of the 
quasar ionizing flux. 



2.4 Formulating the ASC 

Now that we have all the necessary definitions at hand, we 
can describe the formulation of the ASC. As anticipated 
above, the ASC is based on the comparison between i?Hii 
and RT"- 

The first case, Rmi > Rw"^" , corresponds to the PR sce- 
nario, i.e. a situation in which the the quasar photoionization 
rate is too low to keep the gas sufficiently transparent within 
its own H II region. Under these conditions, the measured 
size of the transmission window will be equal to its theoreti- 
cally determined value, R^ = i?™"^, within uncertainties. It 
has already been shown that this situation is likely to occur 
in an ionized universe where quasar H ii regions can expand 
several tens of Mpc into the IGM (Bolton & Haehneh 2007; 
M07), albeit it could also be found in an almost neutral gas 
provided that the quasar luminosity and lifetime are large 
enough to power a rapid growth of i?Hii. 

If instead Rmi < i?™""^ (HR regime), the extent of the 
transmission window is set by the size of the H ii region, i. e. 
by the sharp xm discontinuity at the location of the ioniza- 
tion front. In this case R^^ = i?Hii: by construction, i?™"^ 
represents the upper limit to R^ obtained by assuming that 
quasar ionizing flux extends to infinity. The HR regime usu- 
ally occurs when xm is large {i.e. small i?Hii) and the quasar 
ionizing flux is suppressed abruptly while still dominant over 
the Ultra Violet Background (UVB). If the HR regime holds, 
the measured value R^ gives an estimate of a^ni via eq. [T]as 
first proposed by Wyithe & Loeb (2004). It is worth stress- 
ing that, by their own operational definition, the PR and 
HR regimes are not associated to a particular phase of the 
reionization process. 

Note that, according to our definition, it is always 
i?™"^ ^ Rm, within errors. However, i?™"^ is determined 
under the assumption that Fg = 0, i.e. neglecting a pos- 
sible contribution to F by sources other than the quasar 
and the UVB. This might well not be the case. However, it 
is easy to isolate the effects of contributing sources (most 
likely galaxies) along the LOS which would be signaled by 
a larger than expected value Rw > R^" . In this case one 
would conclude that Fg 0; in addition, it will be also pos- 
sible to obtain a quantitative estimate of such contribution 
by using eqs. |3HlS]and further imposing that F{Rw) = 0.1. 

The application of the ASC will allow to decide between 
the following possibilities: if Rw ^ R^^ the PR regime 
applies and no constraints can be put on xhi; if moreover 
Rw > R^^ this implies a non-vanishing contribution of 
galaxies to the local photoionization rate, Fg 0. If in- 
stead Rw < i?™"^, the HR regime applies, 7?hii ~ Rw, and 
Xm can be constrained. 



3 THE ERROR BUDGET 

In this Section we discuss and quantify the various uncer- 
tainties involved in the determination of the three character- 
istic scales entering the ASC, together with their associated 
error sources. 
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3.1 Rmi 

Rmi denotes the average extent of the quasar H ii region 
and consequently the uncertainty on this quantity is associ- 
ated to fluctuations in the density field and radiative trans- 
fer effects [e.g. shadowing, filtering and self-shielding) which 
distort the shape of the quasar H ii region; the extent of the 
H II region is then direction-dependent. The amplitude of 
this effect can be estimated only with the support of 3D 
numerical simulations. 

We have then performed a set of dedicated 3D ra- 
diative transfer simulations which follow the expansion of 
the If II region in a cosmological density field for differ- 
ent quasars and IGM parameters, as described in more de- 
tails in Sec. 4. We use the outputs of such simulations to 
quantify the relative error Aiinii/iiHii. For each simulation 
we have drawn fOOO LOS towards the quasar and on each 
LOS we have estimated -Rhii as the distance R at which 
dxm/dR > fOO. Ai?Hii is quantified as the 1-a dispersion 
of the -Rhii distribution derived as above. 

Tab. f reports -Rhii and the associated error for var- 
ious simulation runfl We find the relative error on -Rhii 
to be confined in the extremely narrow range 0.07 — O.f f , 
despite the wide range covered by the average extents of 
the H II regions which result from spanning over the values 
allowed for the simulation parameters {N-y, xm, tq). 



3.2 Ru] 

The extent of the transmission window is an observable 
quantity in the spectra of high-2 quasars whose estimate 
is derived from the measurement of the quasar redshift {zq) 
and the onset of the GP trough [z^) along a given LOS. 
The accuracy of the quasar redshift depends on the line in- 
dicators available. The most accurate determination relies 
on CO rotational transition lines from the host galaxy. Only 
for two of the high-z quasars analyzed (see Sec. 5) CO de- 
tection is currently available: Jff48-|-5251 (Bertoldi et al. 
2003) and J0927+200f (Carilh et al. 2007); for these cases 
redshifts are measured with an uncertainty Azq < O.OOOI. 
The next accurate indicator for systemic quasars redshift 
is the Mg II line which, despite being a broad line, allows 
accuracy of the order of 0.002 < Azg < 0.004 (Richards 
et al. 2002). Ifowever the Mg ii transition line has a rest- 
frame wavelength of 2800 A, so that at z ~ 6 it is shifted to 
roughly 2 /^m, thus requiring infrared spectroscopy which is 
not always available. When neither CO and Mg ii detections 
are available the redshift is measured via the Lya line which 
has an associated error of the order of Azq « 0.02. In the 
quasars list reported in Table 1, the redshifts are tagged 
with a letter specifying the line indicator used. 

Extracting an accurate measurement of z^ from the ob- 
served spectra is also not trivial, as already discussed in 
Sec. 2.1. As mentioned there, we adopt the operational def- 
inition given in F06, i.e. z^ is taken as the redshift corre- 
sponding to the first pixel, moving away from the quasar 
location, at which the flux smoothed by a top-hat fllter of 



For runs R2 and R4 these values are not given as these cases 
correspond to an initially highly ionized IGM in which the quasar 
H II region becomes larger than the computational volume. 



20 A drops below the threshold F = 0.1. Atz^G, 20A cor- 
respond roughly to IS.z^ ~ 0.02. We take the latter as a 
measure of the error on z^. 

Aside from the intrinsic error associated to the -R™ mea- 
surement along the observed LOS, the uncertainty on -R™ 
must take into account the dispersion along different lines 
of sight toward the same target quasar. Analogously to the 
case of -Rhii, quantifying such an error requires the support 
of 3D numerical simulations. We have used the same sam- 
ple of LOS used for the A-Rhii determination, to extract the 
mock spectrum and the associated value of -R™ for each LOS 
(the details of such procedure are given in Sec. 4). The av- 
erage -R„ values found for the various simulations are given 
in Tab. 1 together with the associated \-a dispersion of the 
sample. As already pointed out in our previous work (M07), 
-Ru, fluctuations along different LOS are signiflcantly larger 
with respect to the -Rhii ones: this is due to the fact that 
Rw depends on the Lya optical depth which is extremely 
sensitive to tiny fluctuations of xhi, while -Rhii is set by the 
sharp a::Hi transition occurring at the ionization front. Dif- 
ferently from -Rhii, Rw has a stronger dependence on the 
physical state of the gas in which the quasar is embedded. 
We flnd /S.Rw/ Rw values ranging from 0.17 to 0.58, with the 
general trend of AR^/R^ increasing with decreasing values 
of the ratio N^/Tb- This behavior is determined by the fact 
that Rw gets smaller for smaller A^^,, while for higher Vb 
values the dispersion, AT?™, increases as a consequence of 
a more sensitive response to density fluctuations. Given the 
limited a-priori knowledge of N-^/Tb in this study we assume 
the average value {AR^, / Rm) ~ 0.3. 

As described in Sec. 2.3, the determination of -R™"^ involves 
the assumption of the following parameters: the quasar in- 
trinsic luminosity, N-,, the gas temperature, T, and the UVB 
photoionization rate, F-B. 

Quasars luminosities are commonly estimated by con- 
verting the observed magnitudes in number of ionizing pho- 
tons emitted per unit time. Such a procedure requires the as- 
sumption of a specific template for the quasar intrinsic spec- 
trum, which is poorly constrained for high redshift quasars. 
In principle each quasar is characterized by its own spectral 
energy distribution and it would be possible to associate at 
each observed object a power law with a specific spectral 
index inferred from infrared data, which however are not al- 
ways available for the current analysis. In the current work 
we derive A'^ by using the measured absolute magnitude 
-Mi45o at 1450 A, and assuming the spectral template by 
Teller et al. (2002) . This is the composite spectrum of lower 
redshift quasars (z < 3.7), which is well fitted by a double 
power-law with a break close to the Lya line at A = 1280 
A and spectral indexes oeuv = —1.76 and aNuv = —0.69. 
As done in Yu & Lu (2005), we take into account the uncer- 
tainty associated to the template by considering also the 
Telfer double power-low fits to the subsamples of radio- 
quiet and radio-loud quasars, which has spectral indexes 
{aEuv = -1-96, a^riuv = "0.67) and (a^^y^ = -1.57, 
ct^jjv ~ -0.72) respectively. 

The IGM temperature as well as the UVB photoioniza- 
tion rate, Tb, are also poorly constrained at high redshift. 
Nevertheless the statistical analysis of the Lya forest de- 
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tected in high-2 quasar spectra, as well as both analytic and 
numerical theoretical models of the reionization process, al- 
low to bracket the above quantities into a relatively narrow 
range of possible values {e.g. Choudhury & Ferrara 2006; 
Gallerani et al. 2008a). In the application of the ASC to the 
observed high- 2; quasar spectra presented in this paper and 
discussed in detail in Sec. 5, we assume Fa to evolve accord- 
ing to the Early Reionization Model (ERM) developed by 
Gallerani et al. (2008a). In such model, reionization com- 
pletes at Zrei ~ 7 after an extended partial ionization phase, 
necessary to recover the recent determination of the electron 
scattering optical depth found from WMAP5 data (Dunk- 
ley et al. 2009). Such model also simultaneously accounts 
for most of the available data including Lya/Ly/3 GP opac- 
ity, Lyman Limit Systems, cosmic star formation history, the 
number density of high-redshift sources, and the Luminosity 
Functions of Lyman Alpha Emitters at 4.5 < 2 < 6.6 (Dayal 
et al. 2008). Moreover, the UVB photoionization rate pre- 
dicted by the ERM is consistent with the results by Bolton 
et al. (2005, 2007) in the range 4 < 2 < 6. At 2 « 6 the latter 
authors provide a solid upper limit Fs ^ 0.34 x lO"^'^ s~^, 
which we take here to quantify the uncertainties on Fs(2). 
Smaller values of F_b do not affect significantly the R^""^ 
estimate, due to the fact that at the characteristic distances 
involved (several physical Mpc) the local photoionization 
rate is still dominated by the quasar ionizing flux. 

The ERM prescribes IGM temperatures in the range 
T = (1 ± 0.5) X 10* K. We take T = (1.5 ± 0.5) x lO" K, to 
allow for a possible local boost induced by the quasar. Note 
in fact that, as seen from our RT simulations presented in 
Fig. 1, the extra heating provided by the quasar can rise 
the gas temperature above 2 x 10* K only for the case in 
which the quasar ionization front expands in a completely 
neutral gas, a possibility already excluded by various high-2 
experiments (see Introduction). 



4 TESTING THE ASC WITH SIMULATIONS 

Before applying the ASC to observations we have tested it 
through a set of numerical simulations which combine SPH 
(Springel & Hernquist 2003) and full 3D radiative transfer 
(RT) calculations performed with CRASH (Maselli, Ferrara 
& Ciardi 2003). These tests are meant to check the validity of 
the ASC and more precisely to verify that the semi-analytic 
method devised to determine R^""^ gives reliable results. If 
so, when applied to the synthetic spectra extracted from the 
simulations, the ASC should be able to correctly predict the 
physical regime (PR or HR) prevailing in the computational 
volume. 

The simulations setting is similar to the one described 
in M07, to which we defer the reader for details. The two 
main differences here are the inclusion of helium physics 
and the newly developed "colored packets" RT algorithm 
(Maselli, Ciardi & Kanekar 2008) which allows to more re- 
liably compute the effects of hard spectrum sources. This 
allows to extract mock spectra with improved precision. Us- 
ing the SPH density field at 2 = 6.093 as an input, we have 
run CRASH for different combinations of the free parame- 
ters (^^7, tq, xhi). The set of simulations performed is listed 
in Table 1, where we give the free parameters adopted for 
each run, identified by the name in first column, along with 



values of 7?™°^, Rw, their associated l-cr errors and the cor- 
responding regime determined from the simulations. From 
each simulation we build synthetic Lya absorption spec- 
tra along 1000 LOS through the box, smooth the spectra 
at 20 Aas described in M07, and derive the average value 
of J?u). The corresponding is evaluated following the 

procedure described in Sec 2.3 adopting the values of Fa 
and of the specific RT run, as well as the numerically 
computed/LOS-averaged neutral hydrogen density and tem- 
perature profiles. Due to the inhomogeneities of the den- 
sity field and to RT effects, the quasar H 11 regions found 
in our numerical simulations are not spherical, with signifi- 
cant scatter of the H 11 region size along different directions 
{e.g. see Fig. [l}. We then compute -Rhii by averaging the 
value over all the LOS. For all the simulated configurations 
in which the H 11 region is confined inside the computational 
volume (Rl, R3, R5, R6), we find that the average values of 
-Rhii are consistent with the approximated values given by 
eq.[T] within few percent. This result shows photoionization 
equilibrium, assumed in the calculation of _R™"^ (Sec. 2.3), 
to be an excellent approximation. 

Two examples of PR (left) and HR (right) regimes are 
shown in Fig. [T] taken from the runs R3 and R5, respec- 
tively. The maps show H i number density cuts in a region 
of 50 h"^ Mpc (comoving) on a side with the quasar located 
in the right-bottom corner. Overplotted in the maps is the 
averaged temperature profile within the same region. In the 
corresponding bottom panels of the same Figure, the black 
(red) solid lines show the full-resolution (20 A smoothed) 
absorption spectrum along a random LOS to the quasar, 
from which we determine Rw according to the operational 
definition given in Sec. 2.3. The mean transmitted flux de- 
rived from eq. |4] is shown by the solid blue lines, while the 
dashed curves represent the mean transmitted flux for the 
case of an homogeneous density fleld at the mean density, 
i.e. P{A) = 5{1). 

By comparing the analytical transmitted fluxes curves 
to the mock spectra, it is evident that to correctly compute 
the IGM mean opacity at high-2, the proper PDF of the 
density field, and particularly of underdense regions, must 
be taken into account. At high redshift, in fact, r(A = 1) 
grossly overestimates the effective opacity. This occurs be- 
cause fluctuations above a certain overdensity threshold, 
Atfi, start to become opaque due to the progressive increase 
of the average cosmic density with z, which prevents the gas 
at A ^ Ath to remain sufficiently ionized to be transpar- 
ent. The value of Ath decreases gradually with redshift and 
becomes < 1 at 2 ~ 6. 

Fig. 1 illustrates visually the relation among Rw , Rm 
and -Rhii as a function of the different physical conditions 
of the gas in the two regimes. Note that in these plots the 
mock spectra as well as R^ are computed along a single ran- 
dom LOS; -R™"^ is determined analytically, and -Rhii is the 
average value over all the LOS. This choice is meant to re- 
produce the closest analog to the actual experiment. In the 
PR regime (left panel), we find -R™"^ ~ -R^ as expected, 
implying, according to the ASC, that the quasar H 11 region 
should be larger, -Rhii > -R™ ~ Rw . This prediction is con- 
firmed, as seen from the position of -Rhii shown in the panel. 
On the contrary, in the situation in which 7?™"^ > R^ , as in 
the right panel corresponding to the HR regime, Rw truly 
measures the H 11 region radius. Again the simulation results 
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Figure 1. Simulated neutral hydrogen number density maps (upper panels, color bars give lognni values) and corresponding mock 
spectra (lower) along a typical LOS through the quasar for the PR (left, run R3) regime and HR (right, run R5). The maps are 10.1 
physical Mpc on a side with the quasar located in the right-bottom corner. Overplotted are the corresponding mean temperature profiles. 
The lower panels show the full-resolution spectrum (black solid line), the same smoothed at 20 A (red solid), the mean transmitted flux 
derived from eq.|2] (solid blue) and using F = e~'^(^~^) (dashed blue). The values of (measured along the plotted spectrum) and of 
-R™"^ and -Rhii (averaged over 1000 LOS), are shown by the vertical lines, along with the associated l-cr error (shaded band). 



Table 1. Tests of the ASC with simulations 



Run 


Nj,57 








rtmax 


Rw 




Regime 








[10-i3s-i] 


[10^ yr] 


[Mpc] 


[Mpc] 


[Mpc] 




Rl 


0.56 


0.1 


1.25 X 10-^ 


1 


2.51(2.96) 


2.74 ±0.96 


5.32 ±0.61 


PR 


R2 


0.56 


10-4-2 


0.8 


1 


2.30(2.43) 


2.74 ± 1.60 




PR 


R3 


2 


0.1 


1.25 X 10-=* 


1 


4.83(5.35) 


5.50 ± 1.47 


8.51 ±0.80 


PR 


R4 


2 


10-4.2 


0.8 


1 


4.51(4.93) 


5.66 ± 1.86 




PR 


R5 


2 


1 


0.0 


1 


3.4(4.7) 


3.70 ±0.64 


4.26 ±0.31 


HR 


R6 


2 


1 


0.0 


2 


4.2(5.7) 


4.55 ± 0.89 


5.39 ±0.42 


HR 



confirm the validity of the ASC prediction. In this case by 
inverting eq. [T] we obtain a (nominal) value xm= 0.89 ±0.2, 
consistent with the neutral hydrogen fraction, a;Hi= 1, im- 
posed as initial condition in this ruijfj. 

More systematic tests of the ASC have been envisaged 
with the set of simulations reported in Tab. 1, which span 
the allowed values of the free parameters of the problem (A'^, 
tg, xui)- For each run we give here the values of R^"'^ , Rw, 
and their associated l-cr error (see discussion in Sec. 3). We 



^ Rclativistic effects related to finite light propagation on the 
apparent shape of the H II regions do not affect the direction 
along the LOS and thus the validity of the ASC which applies to 
single LOS. 



provide two different estimates for R^°'^: the first (fiducial) 
is obtained by using the average nni profiles extracted from 
the RT simulations, the second, R^"'^ (noRT) (in brackets) 
by neglecting RT effects. In this latter case xm inside the 
H 11 region is calculated by assuming that the SPH simu- 
lated density field is at photoionization equilibrium with the 
unattenuated quasar flux plus background: this method does 
not account for (i) residual H i inside the H ii region, and 
(ii) the drop of the quasar fiux beyond the ionization front. 

From the data in Tab. 1 we find that J?™"^ is within the 
range derived from the mock spectra; R^"'^ (noRT) is in 
agreement with i?™ for PR cases, but it exceeds Rw for the 
R5 and R6 (HR regime) runs, as expected from the ASC: 
Rw in these cases corresponds to the H ii region radius (eq. 
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[T]), a result confirmed by the analysis of these simulation 
outputs. 

The comparison between the fiducial and no-RT (brack- 
eted) i?™"^ values shows furthermore that RT effects asso- 
ciated to the internal residual H i opacity can be safely ne- 
glected in the calculation of R^"'^, whereas the flux drop 
occurring in the HR regime beyond the ionization front can- 
not be disregarded. Stated in a different manner, the power 
of the ASC to discriminate between PR and HR regimes ex- 
ploits the strong flux drop-off taking place at the H ii region 
ionization front. 



5 APPLICATION TO DATA 

Having successfully tested the ASC, we apply it to a sample 
of high redshift quasar spectra taken from F06. The selected 
quasars name, redshift, ionizing photon rate, observed 
value along with our flducial i?™"^ estimate, are reported 
in Table 2. The redshift values are taken from FOG, un- 
less measurements of low-ionization lines like Mg ii (Kurk 
et al. 2007; Jiang et al. 2007) or of CO rotational lines 
(J0927-h2001 from Carilh et al. 2007; J1148+5251 from 
Bertoldi et al. 2003) do exist. The QSO ionizing rates have 
been derived from their magnitudes, as discussed in Sec. 3.3 
assuming the Telfer composite spectrum. The values 
are determined from the semianalytic procedure described in 
Sec. 3.3 and adopting the overdensity PDF extracted from 
the SPH simulation described in Sec. 4. In addition, Fs 
and the IGM temperature, T, are assumed to evolve ac- 
cording to the Early Reionization Model (ERM) developed 
by Gallerani et al. (2008a). The details of the above mod- 
els have been given in Sec. 3.3, where we also discuss the 
uncertainties associated to Fs and T. 

Fig. [2] shows the comparison between R^^ and R^"'^ 
for the quasar sample in Tab. 2. As discussed in Sec. 3.2, 
the errors on R^ are dominated by the uncertainty associ- 
ated to the dispersion along different LOS. The 7?™°"" val- 
ues shown in Fig. 2 by the black squares are obtained by 
adopting an IGM temperature T = 1.5 x 10*K. The shaded 
region represents the error on f?™"^ associated to the uncer- 
tainties on the temperature and on the spectral template: 
the upper (lower) solid line connects the values associated 
with the larger (smaller) photoionization rates correspond- 
ing to the double power-law fltting the radio-quiet (radio- 
loud) composite spectrum, and to the upper (lower) limit 
on the temperature consistent with the ERM, T = 2 x 10^ 
K {T — 10^ K). Note that higher temperatures at fixed 
density yield lower optical depths, larger transmitted fiuxes 
and consequently higher i?™"^ values. Fig. 2 also shows the 
j^max yg^jygg obtained for the fiducial T and A'^^ values, but 
with different assumptions on Fs: we consider the extreme 
cases of (i) a constant Tb equal to the upper limit given by 
BH07 of Fs,_i2 = 0.34 (dashed-black line), and (n) Fs = 
(dotted- black line). The small differences among these two 
extreme cases demonstrates that uncertainties in the inten- 
sity of the UVB do not sensibly affect the ASC. This results 
from the fact that, at the scales relevant to ASC, the local 
ionization rate is dominated by the quasar flux, i. e. F ~ Fq. 

Fig.[2]represents the central result of this paper, i.e. the 
application of the ASC to a sample containing the most dis- 
tant quasars, and it allows to draw a number of interesting 



conclusions. None of the analyzed quasars satisfles the con- 
dition Rw < (within uncertainties), indicating that all 
sample QSO are in the PR regime. As a consequence, cur- 
rent observations cannot provide reliable estimates of the 
IGM ionization state by studying their H ii region. 

Secondly, we find Rw > i?™"^ for the quasar 
J1602-I-4228 at > 2 — cr significance level, and for the quasars 
J1436-h5007, J1411+1217, J1623-f3112 at a lower signifi- 
cance level. According to the ASC this occurrence signals 
a locally enhanced photoionization rate. As discussed in 
Sec. 2.4 this could be due to either a higher-than-average 
number density of galaxies clustered around the quasar or 
a local enhancement of more distant galaxies along the ob- 
served LOS, but not directly associated to the overdensity 
in which the quasar is located. Another possibility is that 
the contribution to the local photoionization rate due to 
the galaxies around the quasar is increased by the enhanced 
mean free path produced by the quasar ionizing fiux. In 
order to decide among these possibilities accurate spectro- 
scopic redshift determinations from deep galaxy surveys of 
the quasars fields would be necessary, which unfortunately 
are not available yet. 

Independently of its origin, the ASC allows to directly 
infer the amplitude of Fg at R^ by progressively increas- 
ing F(r) (eq. O untfl the condition i?™""" = R^ is met. 
Following this procedure we obtain: Tq'^^^ = 1.63 x 10"" 
s-\ F;^"" = 1.85 X 10-13 s-\ TU^''^ = 2.38 x lO'^^ s'^ 
and Vq^"^^ = 3.09 x 10~" s'V These values are compara- 
ble to the average intensity of the UVB at the same redshift: 
0.35 < Tg/Vb < 0.85. Furthermore, for aU these four QSOs, 
we find that Fg is comparable or exceeds Fq at R^^, with 
derived values of the ratio Fg/Fq = 1.25,1.05,1.49,1.02. 
Such ionizing photon contribution from galaxies is crucial 
to explain the width of the transmission windows. 



6 DISCUSSION 

We have proposed the Apparent Shrinking Criterion (ASC), 
a simple but yet robust semi-analytic method, based on the 
apparent shrinking effect (Maselli et al. 2007), which can 
be applied straightforwardly to the interpretation of trans- 
mission windows in the the spectra of distant quasars. The 
ASC is based on the comparison between three spatial scales 
associated with a quasar: these are Rmi, the average physi- 
cal radius of the H ii region, the extent of the transmission 
window detected in its spectrum, R^n, and R^°'^ , the dis- 
tance from the quasar at which the mean transmitted ffux 
drops to = 0.1 in the absence of attenuation due to either 
H I internal to the H ii region or the ionization front of the 
latter. 

The ASC applied to the spectra of single high-2: QSOs 
provides a novel method to discriminate among the two sce- 
narios previously proposed in the literature to explain the 
QSO transmitted fiux. In the first one the extent of the ob- 
served transmission window, R^ , is set by the sharp drop of 
the quasar ionizing fiux at the ionization front location, so 
that Rw gives a measure of Rmi- If the quasar H 11 region 
grows large though, the unattenuated/geometrically diluted 
quasar fiux, gets progressively weaker and eventually falls 
below the threshold necessary to keep the IGM transpar- 
ent to Lyof radiation. As a result, the observed R^ satu- 
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4.40 


J0002 


+ 2550 


5.80'' 


5.650'' 


—27.70 


2.33 


1.98 


2.68 


10.68 


7.15 


J0836 


+ 0054 


5.81'= 


5.721'* 


-27.83 


2.63 


3.02 


2.23 


6.24 


7.55 


J1436 


+ 5007 


5.83'' 


5.720'' 


-26.51 


0.79 


0.90 


0.66 


7.69 


3.75 


J0005 


- 0006 


5.85^* 


5.771'* 


-26.42 


0.72 


0.82 


0.61 


5.45 


3.55 


J1411 


+ 1217 


5.93^= 


5.823'* 


-26.70 


0.93 


1.07 


0.80 


7.20 


3.85 


J0818 


+ 1722 


6.00'' 


5.920'' 


-27.37 


1.72 


1.98 


1.46 


5.22 


4.95 


J1306 


+ 0356 


e.oi'* 


5.916'* 


-27.14 


1.39 


1.60 


1.18 


6.13 


4.45 


J1137 


■f 3549 


6.01'' 


5.910'' 


-27.08 


1.32 


1.51 


1.12 


6.53 


4.25 


J1630 


+ 4012 


6.05'' 


5.980'' 


-26.11 


0.54 


0.62 


0.46 


4.48 


2.65 


J1602 


+ 4228 


6.07'' 


5.950'' 


-26.82 


1.04 


1.19 


0.88 


7.70 


3.65 


J1250 


+ 3130 


6.13'' 


6.030'' 


-27.11 


1.35 


1.56 


1.15 


6.25 


4.05 


J1623 


+ 3112 


6.25^= 


6.160^ 


-26.67 


0.90 


1.04 


0.77 


5.39 


3.05 


J1030 


+ 0524 


6. SI'* 


6.217'* 


-27.15 


1.41 


1.62 


1.20 


5.52 


3.65 


J1148 


+ 5251 


6.42<= 


6.325^ 


-27.82 


2.60 


2.99 


2.22 


5.37 


4.75 



Table 2. Quasars Sample. Quasar redshifts are tagged according to the spectral lines used for their measurement and to reference works 
as follows: (a) CO (Carilli et al. 2007), (b) Lya (Fan et al. 2006), (c) Mg II (Jiang et al. 2007), (d) Mg II (Kurk et al. 2007), (e) 
CO (Bertoldi et al. 2003), (f) (White et al. 2003.) 




Quasar redshift 



Figure 2. Comparison between the observed and corresponding values of i?™"^ for the sample quasars listed in Tab. 2. Red 
diamonds show the measured R^, black squares are the values obtained for ij™"^ . The shaded region accounts for the errors associated 
to the uncertainties in the gas temperature and quasar template; the dashed (dotted) line connects i?™"^ values obtained assuming 
Tb, -12 = 0.34 (Tb = 0), at fixed IGM temperature (T = 1.5 X 10*) and Teller quasar template. 
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rates at the theoretically estimated upper limit R^°'^ , in- 
dependently of the actual extent of the H ii region, so that 
Rw ~ -Rm"^ < -Rhii. We refer to the first/second case as 
the HII region (HR) /Proximity Region (PR) regime, respec- 
tively. In the HR regime, one can determine xm, the IGM 
mean H i fraction. If the PR regime holds instead, the value 
of Rw can be used to constrain the local (i.e. at a distance Rw 
from the quasar) intensity of the ionizing radiation and indi- 
rectly infer the contribution of intervening clustered galaxies 
to the local photoionization rate. 

As formulated in this paper, the ASC can only be ap- 
plied to single objects and at the moment it is unclear 
whether it could be extended to get information on the av- 
erage IGM properties at high redshift. At the moment, the 
ASC must be interpreted as a useful local test of reioniza- 
tion, i.e. in regions surrounding luminous quasars. Neverthe- 
less, by separately collecting information from a large sample 
of quasars in a given redshift range, the ASC can eventually 
give new constraints on the global reionization history. We 
leave the analysis of this further extension of the ASC for 
future work. 

The novelty introduced by the ASC is that it makes pos- 
sible to assess the physical regime (HR or PR) characterizing 
an observed quasar. Previous works have instead focused on 
extracting global information on the IGM redshift evolution 
from generic trends in the data, e.g. the dependence of Rw 
on redshift (FOG; W08), statistical fits to the Rw distribution 
function (Wyithc & Locb 2004), and/or fits to the flux pro- 
files within the transmission windows (Mesinger & Haiman 
2007). The drawback of the first two approaches is that, in 
order to carry such studies, it is necessary to make a scaling 
of Rw to a reference quasar luminosity, a procedure which 
is valid only in the case in which the quasars are all in the 
HR regime. By implicitly making this assumption the infor- 
mation encoded in the observed spectra is degraded. The 
same happens in the flux profile fitting procedure, which is 
meant to constrain the most probable values for the triplet 
(-Rhii, N.y, xhi). Although in this case the hypothesis that 
Rw identifies the H ii region radius is not made explicitly, 
the transmission profile shape is sensitive to the three pa- 
rameters above only in the HR regime; in the PR regime the 
profiles arc sensitive to A^-y, but completely degenerate with 
Run, and xm. 

The Apparent Shrinking Criterion has been extensively 
and successfully tested, by applying it to synthetic absorp- 
tion spectra extracted from a set of state-of-the-art numeri- 
cal simulations. Next the ASC has been applied to a sample 
of high redshift quasar spectra taken from Fan et al. (2006). 

The main result is that all quasars in the analyzed sam- 
ple arc found to be in the PR regime and hence their spectral 
observables (inner flux profile, extent of transmission win- 
dows) cannot be used to reliably constrain the value of the 
average IGM neutral hydrogen fraction. Nevertheless, for 
quasars in the PR regime the ASC offers a novel method to 
determine the local photoionization rate, analogous to the 
proximity effect at lower redshift (Bajtlik et al. 1988). As an 
outcome of our analysis, we find that all the observed quasar 
spectra considered are consistent with the Early Reioniza- 
tion Model described in Gallcrani ct al. (2008a), which pre- 
dicts reionization to be complete at a « 7 and thus the IGM 
at ^: w 6 to be highly ionized. 

Beside these general results, we have found that 



four objects (J1436+5007, J1411-hl217, J1602+4228 and 
J1623+3112) show the presence of an additional contribu- 
tion. To =^ 0, to the local photoionization rate. The origin of 
this extra ionizing radiation is still unclear. Such flux could 
be provided by (i) galaxies clustered in the biased high den- 
sity regions in which quasars are expected to form, (ii) a 
galaxy overdensity intervening along the observed LOS to 
the quasar and not directly associated to the quasar over- 
density, (iii) galaxies within the quasar H ii region whose 
contribution to the photoionization rate is boosted by the 
quasar-enhanced mean free path of ionizing photons. For 
the four quasars above we found that 1.02 ^ Fg/Fq ^ 1.49, 
thus showing that for these cases the galaxy ionizing flux 
is comparable to the quasar one at Rw- The interpreta- 
tion of this result within structure formation scenarios might 
give interesting clues on the environmental impact of high-z 
quasars. 

In the process of testing the ASC with our RT-I-SPH 
simulations, we have shown a further interesting point: a 
correct estimate of the IGM mean opacity at high- 2; requires 
to properly weight the contribution from underdense regions 
with their actual probability distribution function. At high 
redshift in fact, fluctuations with overdensity above a certain 
threshold start to become opaque due to the progressive 
increase of the average cosmic density with z, and hence do 
not contribute to transmission even if they are still highly 
ionized. This overdensity threshold decreases gradually with 
redshift and it falls significantly below unity at 2 ~ 6. 

As a final remark, we emphasize that the lack of sam- 
ple quasars in the H 11 region regime, does not dismiss the 
potential use of quasar transmission windows to directly 
measure xhi at 2 > 6: in fact, we do expect that with the 
detection of higher redshift and/or lower luminous objects 
expected in the near future, spectra in the HR regime will 
likely be found. 
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